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Abstract 

The processes 77 — > vr^vr'^Tr'^ and 77 tt^tt^tt^ are considered in Generalized Chiral 
Perturbation theory, in view of their potential sensitivity to the mechanism of spontaneous 
breaking of chiral symmetry and to various counterterms. The amplitudes are computed 
up to order 0{p^). The event production rates are estimated for the Daphne (/>-Factory 
and for a future r-Charm Factory. 
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1 Introduction 



In the limit where the masses of the hghtest quark flavours u, d and s are set to zero, the 
QCD lagrangian becomes invariant under a chiral SU{3)i^ x SU{3)b. global symmetry. This 
symmetry is not reproduced by the hadronic spectrum, and must therefore be spontaneously 
broken towards the diagonal SU{3)y subgroup. Actually, this spontaneous breaking of chiral 
symmetry can be shown to follow from very general properties of QCD [|T], 0, |^. However, 
besides the existence of eight massless pseudoscalar states coupling to the eight conserved axial 
currents, nothing is known from "first principles" about the actual mechanism of spontaneous 
breaking of chiral symmetry. The widespread belief in that matter is that it proceeds through 
the formation of a strong quark-antiquark condensate, < qq >~ — (250Mey)^, where < qq > 
denotes the single flavour quark-antiquark condensate in the SU{2) chiral limit, 

< qq> = <uu> |^^=^^=o = <dd> |^^=^^=o ■ (1-1) 



In particular, this means that once the light quark masses m^, and are turned on, the 
mass of the pion is assumed to be dominated by the contribution linear in quark masses 0], 

^ , m = , (1.2) 



where F stands for the pion decay constant (the normalization we use corresponds to the 
numerical value = 92.4 MeV [^) in the same two-flavour chiral limit, 

F = F7r|m„=m^=0- (1-3) 



However, our present theoretical knowledge of non-perturbative aspects in QCD does not 
exclude a picture where the condensate would be much smaller, say < qq >~ — (lOOMeK)^, or 
even vanishing. How the latter possibility may arise in QCD has been discussed recently 
in terms of spectral properties of the Dirac operator. On the other hand, it has also been 
suggested [§] that a strictly vanishing condensate in the chiral limit could be excluded by an 
inequality between the correlator < Afj_{x)A'^{0) > of two axial currents and the two-point 
function < P{x)P{0) > of the pseudoscalar quark bilinear density P{x) = {qil5q){x). This 
inequality, however, has only been established so far for bare quantities, i.e. in the presence 
of an ultraviolet cut-off Ayy. As the cut-off is removed, the pseudoscalar density P{x) needs 
to be renormalized, whereas the (partially) conserved axial current A^{x) remains unaffected. 
Stricktly speaking, the claim of Ref. p is that < qq > (Auv) 7^ for finite Auv, which does not, 
a priori, exclude the possibility of a vanishing condensate in the limit Auv 00 H. Finally, 



Ref! . [10| provides an example of a different approach which, within a variational framework, 



leads to a non-vanishing, but nevertheless small, value of < gg >. 
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Clearly, in order to settle the question of the mechanism of spontaneous chiral sym- 
metry breaking on a purely theoretical level, a major breakthrough in our understanding of 
non-perturbative aspects of confining gauge theories is required. Instead, one may try a phe- 
nomenological approach, and look for experimental observables which could provide the relevant 
information. In this respect, low-energy tt-tt scattering in the S-wave has been put forward as 
a process particularly sensitive to the size of the condensate [|lT], [1^ (for recent discussions, see 



e.g. HTSl , 1^). This can be most conveniently seen in the framework of the effective lagrangian 



T^, |T6|. In order to incorporate the possibility of a small condensate, the usual counting 
has however to be modified. As explained in Refs. Jll], |17|, a consistent expansion scheme, 
usually refered to as Generalized Chiral Perturbation Theory (GChPT), is obtained by taking 
(from now on, we restrict ourselves to the case of two massless flavours only), 

m^,mdr^Oip) , B^O{p) , (1.4) 

with p being a generic momentum, much smaller than the typical hadronic scale Ah ~ 1 GeV, 
and 

B^-^. (1.5) 

With this counting, the effective lagrangian at lowest order, which consists of all chiral invariant 
terms of order O(p^), reads 

+A{{U+X? + ix^UY) + Z^iU+x - X^Uf (1.6) 
+ hoix'^x) + hi{detx + detx^)} . 

The matrix U collects the pion fields (throughout, we adopt the Condon and Shortley phase 
convention). 



The notation is as in Refs. [^, |I6[, except that x^ the quantity that contains the scalar and 



pseudoscalar sources, is defined without the usual factor 2B, 

X = s + ip = M -\ , M = diag(m„, nid) . (1.8) 

The covariant derivative contains the external vector and axial sources, 

D^U = d,U - t[v^, U] - i{a^, V} . (1.9) 
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This is to be contrasted with Standard Chiral Perturbation Theory (SChPT) |1T5| , 16], which 
assumes a large condensate, and hence the counting rule 171^,171^ ~ C^(p^)) B ~ A^, so that 
only the first two terms on the right-hand side of Eq. ( |1.6|) are taken into account at leading 
order. It easily follows from ( |1.6|) that the lowest order expression of the pion mass is now given 

byQ 

= 2mB + Am^A + ■ ■ • , (1.10) 
where the ellipsis stands for higher order corrections. At the same level of approximation in 



the chiral expansion, the tt-tt scattering amplitude can be written as |11 



A{s%u) = p ^' ;r + ' (1-11) 

with P = 1 + 0(jri), while at this order the parameter a is directly related to the < qq > 
condensate through 



a 



The case ( |1.2|) of a strong condensate corresponds to a = 1, whereas the extreme limit where 
the condensate would vanish yields a = 4. Notice that the above expression for A{s\t,u) 
is not affected by the additional terms in ( p..6|) and reproduces the result first obtained by 
Weinberg ]18|], A{s\t,u) = {s — 2mB)/F'^ + ■ ■ ■. The difference between the standard case 
and deviations from it lies here only in the leading-order expression ( |1.1CID of the pion mass 
and its relation to the condensate. Higher orders in the chiral expansion will modify the 
simple expression (|1 . 1 1|) , but the correlation between low-energy tt-tt scattering and the value 
of the ratio 2mB/M^ subsists, and can be studied in a controled way within the generalized 
chiral expansion [0. Available data on low-energy tt-tt phases, which are dominated by the 



Geneva-Saclay experiment ||T^, do however not possess the required accuracy in order to 
distinguish between the different alternatives at present. Forthcoming experiments, such as 
new experiments, conducted by the BNL865 collaboration [?] or planed at the Daphne 



factory |2^, and the DIRAC experiment at CERN [^, represent promissing prospects in 



this direction. 

In order to create a (hopefully convergent) set of evidences pro or contra a specific picture 
of spontaneous breaking of chiral symmetry, it remains however important to explore other 
possibilities, and to find other processes which, at the theoretical level, can be shown to exhibit a 
reasonably strong dependence on the value of the condensate. The present work was motivated 
by the above consideration and the following observation. A straightforward re-analysis in 
GChPT of the existing SChPT calculations |2^, at lowest order shows that the two 



^We neglect the mass difference r7i„ — and set m„ = rrid = fa. 
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amplitudes for the production of three pions in low-energy photon-photon collisions involve 
the parameter a already at tree level. In the case of 77 — > tt'^tt^tt^, the amplitude reads 



A 



N 



47r2F3 



s-Ml 



+ 



while for 77 — > 7r"'"7r 7r° we find 



^2 



47r2F3 



k' ■ p_ 



s-Ml 

p_ 



i:^ P-) + 



:i.i3) 



(1.14) 



In these expressions, A; and k' {e and e') denote the momenta (polarizations) of the two photons, 
while p+, p_ and pq are the pion momenta. The ellipses stand for higher order corrections which 
will be considered later. Diagrammatically, the origin of the dependence on a lies in the presence 
of one-pion reducible contributions to the two amplitudes, see Fig. 1 below. 






Figure 1 : The lowest order contributions to the 77 — > tttttt amplitudes. The a-dependence in 
and A^ comes from the vertex for (virtual) tt^ scattering in the first graph, which is the 
only one to contribute in the case of the neutral amplitude. 



Therefore, the leading-order neutral amplitude A^ increases as the condensate decreases. For 
a strictly vanishing condensate (a = 4), the cross section at low energies is thus enhanced by 
a factor 16 as compared to the standard case of a strong condensate {a — 1) \ In the charged 
case, the situation is less favourable. At threshold, the average over the photon polarizations 
of the modulus squared of the amplitude is still proportional to the square of a, 

but the sensitivity on a of the corresponding total cross section a'"(s,Q;) rapidly decreases 
with increasing energy. For instance, the ratio (T^{s,a = 3)/a^{s,a = 1) is equal to 4.48 at 
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a/s = 450 MeV, i.e. just above threshold, but drops to 1.68 at i/i = 500 MeV and becomes 
less than 1.10 at ^/s ^ 600 MeV. For the case a = 2, the corresponding ratio is equal to 2.21 
at ^/s = 450 MeV, but the effect is less than 25% at = 500 MeV, while it barely reaches 
a few percent at y/s > 600. Thus, the interference between the two kinematical structures 
contributing to the amplitude A'" , which is responsible, at low energies, for the suppression 
of the cross section in the charged channel as compared to the cross section in the neutral 



channel ||2^, also washes out the dependence on the value of the < qq > condensate as soon 
as one leaves the threshold region. Considering definite polarization configurations for the two 
photons does not improve the situation : For parallel polarizations of the two photons, the part 
of A'" which is sensitive to a does not contribute to the cross section, whereas for polarizations 
taken along orthogonal axes, the same interference effect is again fully at work. Furthermore, 
both cross sections rapidly rise above threshold, so that the effect of higher orders also needs to 
be investigated. The purpose of this paper is precisely to investigate these higher order effects 
at the one-loop level in GChPT. In fact, at next-to-leading order new tensorial structures 
appear both in the neut! ral and in the charged amplitude p5[, whereas the pion loops provide 



additional sources of dependence with respect to a, so that their behaviour with respect to 
changes in the value of the condensate could be modified to some extent. 

Accordingly, the rest of the paper is organized as follows : The general kinematical struc- 
ture of the two amplitudes A'^ and A'^, as well as their properties under isospin symmetry, are 
the subject of Section 2. The construction and the renormalization of the effective lagrangian 
of GChPT at order one-loop in the two-flavour case are treated in Section 3. Section 4 is de- 
voted to the actual calculation of the two amplitudes to one-loop precision. The counterterms 
which are involved in these expressions and several numerical results are discussed in Section 
5. A summary and conclusions are presented in Section 6. Details on the evaluation of some 
counterterms have been gathered in an Appendix. 



2 Kinematics and Isospin Symmetry 

The amplitudes A^ and for the processes 

j(k)^ik') n\p,) n\p,) n\p,) , 

^{k)^{k') ^ 7r+(j)+) 7r-(p__) ttV) , (2.1) 
are obtained from the matrix elements 

7r°(pi)^°(P2)vr°(p3) out I 7(A:, e^k', e') in > = z(27r) V(P; - P,)^^ , 
■-(p+)7r-(p_)7r°(po) out I 7(A;, e)^{k\ e') in > = 2(27r) V(P/ - P,)^^ , (2.2) 



< 



respectively, with 



A {k,e;k',e';po,p+,P-) = 

ie\ik)elik') [d''xe~'^-^< 7r+(p+)7r-(p_)7rO(po) owt | T{j^(x)/(0)} | > .(2.3) 



In the above expressions, J^(x) denotes the electromagnetic current, with its usual decompo- 
sition into an isotriplet and an isosinglet component, = + while | Q > stands for the 
QCD vacuum with massive light quarks, but in the absence of electromagnetism (radiative 
corrections to the two processes (|2.1|) are not considered here). Since we also neglect isospin 
breaking effects due to 7^ m^, Bose symmetry and G-parity constrain the three final pions 
to be in an / = 1 total isospin state. Thus the matrix elements in Eq. (|2.3|) only involve the 
1 = 1 component of the product of the two electromagnetic currents. 



T{Ux)UO)}i=i = T{jUx)fM+f,ix)fM}. 



(2.4) 



Therefore, isospin invariance relates the two amplitudes A'" and A^ in a simple way. With the 
Condon and Shortley phase convention adopted in Eq. ( |1.7| ), this relation reads 



A^{k,e; k',e';pi,p2,P3) = A^'{k,e; k' , e';pi,p2,P3) + cyclic {pi,P2,Pz) 



\C 



(2.5) 



where "cyclic (pi,P2;P3)" indicates that the contributions arising from cyclic permutations over 
the pion momenta pi,p2 and pz have to be added. 

Up to permutations of the momenta and polarizations of the photons, and/or permuta- 
tions of the momenta of the charged pions, the amplitude A^ may be decomposed into six 
independent Lorentz invariant amplitudes 



A {k,e;k',e;po,p+,p-) = A^{k,k';po,p+,p^)ti{k,e;k',e') 



+ E 



=2,3 



A ik, k'; po, p+, p_) ti{k, e; k', e'; po, p+, p_) + 



'k^ 



(2.6) 



+ E 

j=4,5,6 



A^{k, k'; Po, p+, P-) U{k, e; k', e'; po, p+, P-) + 



.e' 



+ 



P+ ^ P- 



with {pij =pi+ pj, i,j = +, -, 0) 

ti(A;,e;fc',e') = e^,„^Pe'^A;'V^ 



hik,e; k\e';po,p+,P-) = e,.„^fc"^6>?(-6^ + ^^pi + f^/) 



k ■ p- 



k-p+ 



ts{k,e; k',e';po,p+,P-) = e,,^p{-e^ + l^p) fc^^e'/^. 



k-p+. 
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t5{k,e;k',e;po,p+,P-) 



U{k,e; k',e';po,p+,P-) 




) 



(2,7) 



In the sequel, we shall compute the amplitudes A'" and A^ within the framework of 
generalized chiral perturbation theory up to order one-loop. At this level of accuracy of the 
chiral expansion, the amplitude A^ does not yet receive its full structure as implied by Eqs. 



where the ellipsis stands for higher order terms in the chiral expansion. In order that ( |2.8| ) 
follows from (|2.5|) and ( pl6|) , it is sufficient that the one-loop charged amplitudes A^ satisfy 
the following conditions : 

i) ^^(A;, A;';po)P+;P-) ~ •^5{k,k';p+,p_,po) is entirely symmetric under permutations of the 

pion momenta p+, P-,Pq; 

ii) k ■ po[A'^{k,k';po,p+,P-) + A^{k,k';po,p+,P-)] = 



ill) A^{k,k';po,p+,P-) = A'^'{k,k';pj,,po,P-). 

Furthermore, the two neutral amplitudes have then the following expressions in terms of the 
charged ones : 




(2.8) 



k-p+ [A^{k,k';p+,po,p-) + A^{k,k';p+,po,P-)]; 




{k,k';pi,p2,P3) 




(2.9) 



and 




ik,k';pi,p2,P3) 



[A^{k, k';pi,p2,P3) +AQ{k, k';pi,p2,P3)] + [pi ^ P2] 
+A^ik,k';pi,p2,P3) + A^{k,k']p2,PuP3) ■ 



(2.10) 
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3 The Effective Lagrangian in GChPT 



Before proceeding with the calculation of the amplitudes A*" and A'^ in the next section, we first 
discuss the structure of the low-energy generating functional in the case of two light flavours 
that we need for our subsequent calculation. Since most of the results of this section are not 
available from the existing literature, we discuss them in some detail. 

The structure of the effective lagrangian C^^ is independent of the underlying mechanism 
of spontaneous chiral symmetry breaking. It consists of an infinite tower of chiral invariant 
contributions 

^'"^ = E ^m) , (3.1) 

where C(^k,i) contains k powers of covariant derivatives and / powers of the scalar or pseudoscalar 
sources. In the chiral limit, these terms behave as 

with TTiquark = ^uifT^di and p stauds for a typical external momentum. The standard approach 
not only assumes mquark ^ A//, but also mquark ^ B/2A, such as to enforce the dominance of 



the term linear in m-quark in the expression of the pion mass ( 1.10 ). This allows to reorganize 



the double expansion ( pTLf ) as [|15[ |r 



£cfr ^ ^(2) _^ ^(4) _^ ^(6) ^ . . . ^ (3_3) 



where C^^^ = ^{k,i) with k + 21 = d. The generalized framework considers the possibility 
that the condensate could be much smaller than usually believed, so that for the actual values 
of the quark masses one could have mq^ark ~ B/2A and still mq^ark ^ ^h- This leads to a 
different reorganization of the double expansion ( p.l| ), namely 

£efr ^ ^(2) ^ ^(3) ^ ^(4) ^ ^(5) ^ ^(6) ^ . . . ^ (3_4) 



where now, according to (|r|), C'-'^^ = E B''C^k,i) with k + l + n = d^,^. 

The leading order of the generalized expansion is described by C^'^\ which in the two 
flavour case was given in Section 1, Eq. ( [L.6| ). For our purposes, we need only to consider the 
situation without axial source, and with the vector source restricted to the (classical) photon 
field A., 



"A" 







eA.Q 



(3.5) 
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where Q stands for the charge matrix of the two hght quark flavours u and d, 




(3.6) 



In GChPT, the next-to-leading-order corrections are of order (9(p ), and still occur 
before the loop corrections. They are embodied in C^^^ = C(2,i) + -^(o.s)! which reads ^ 

^(3) ^ lF'[^i^){D,U+D^U{x'U^U+x)) 

+pf (x+f/ - U+x){{x^Uf - [U+xf) (3.7) 

/^From order 0(p'^) onward, the contributions to C^^ come with either even or odd intrinsic 
parity, jC^'^'> = for c/ > 4, or C(^k,i) = '^(fe i)+^lk i) > 4. The tree-level contributions 

at order C(p^) in the even intrinsic parity sector are contained in 

The part without explicit chiral symmetry breaking, £^ q-j , is described by the same low-energy 



constants /i, /2, ^5, h and /12 as in the standard case 15 



- (2/^2 + Ih) ( F^^.F^^'^ + F^.F^'^^ ) . (3.9) 
The part with two powers of momenta and two powers of quark masses is given by 

Ci2,2) = ^F^I^m{D,U+D^^Uix^x + U+XX^U)) 
+a2{D^U+Ux^D^UU+x) 

+a,{D^U+U{x^D^^X - D^'X^X) + D.UU^ixD^^x^ - D'^XX^)) 

+h{D,U+D^^Uix^Ux^U + U+xU^x)) 

+b2{D,U^xD^ U^X + X^D^Ux'-D^U) 

^ The superscript (2) is meant to distinguish the low energy constants ^^^^ , p^"* from the similar ones that 
occur in the expression of C'^^^ in the three flavour case lO, O. 
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+ci {D^U^x + X'D^U) {D^^U^x + X+^^f/) 
+C2(D^X+f/ + U+D^x){D^^U+x + 

+C3(D^X+t^ + U+D^x){D^^X^U + fZ+Z^'^x) (3.10) 
+C4(D^f/+X - X+/^;.f/)(/^^f/+X - X-'D^U) 
+c,{D,x^U - U^D^x){D^X^U - U^D^x) 

Finally, the tree-level contributions which behave as C'(^quark) chiral limit are contained 

in >C(o,4), which reads 

+e2(x+x(x+f/x+f/ + f/+xf/+x)) 
+ez{x^xU^XX^U) 

+f2{ix^u)' + (t/+x)')(x+f/ + u+x) 
+/3(x+x(x+f/ + f/+x)) (x+f/ + u+x) 
+U{{x'uf + (f/+x)')(x+t/ + u+x? 
+Mix^uf - iu+xf){x^u - u^x) 

+hA{x'^XX'^x) 
+/^5(x^x)(detx + detx+) 

+/i6(detx + detx^)^ 

+/i7(detx-detx+)'} • 

Notice that in the standard framework the contributions from >C{2,2) and from £(0,4) would count 
as order O(p^) and order (9(p^), respectively 0. 

Next, we turn to the odd intrinsic parity sector. There, the first contribution starts at 
order (9(p^), and since it is entirely fixed by the short distance properties of QCD, there is no 
difference between the standard and the generalized case, = = In the two flavour 

case, £^ vanishes in the absence of external sources. In the presence of an electromagnetic 
field, it reads 

/:^4,o) = ^^e^''''^A,{Q{d,udo,u^df,uu+ - d,u+d^udpU+u)) (3.11) 

,>2 1 1 

-—^e^'''''^d^A,A^{Q^dpUU+ + Q^U+d^U - -QUQdpU+ + -QU+Qd^U) . 



' The standard C'(p^) effective lagrangian = £(6,0) +'^(4,1) + ^{2,2) +^^(0,3) has been worked out in Ref. 
for the case of three light flavours, and very recently, for both two and three light flavours, in Ref. p8|. 



10 



The computation of the amphtudes and A'" to one loop also involves the counterterms 
from C^^^ = and from = C^qq^ + ^^^42)- -'■^ standard case, both '^^(^ 1) and C'^qq^ 

count as order 0{p^), and have been discussed before in the literature I^O], |31|, Borrowing 
from the last and most recent of these references, we obtain 



(4,1) - 4^2 

+Au{[D^UUD''UU[x]-[G''^ 



+Ars{[D'^UUx]-[D''UUG''^U) 



(3.12) 



and 



+iA3{[D''U]4[D^G""']4G''^]+ - 

+A8([/^"/^^f/]-([D'^f/]_[D^[/]_[G'n+ - [G'n+[^7f/]-[^^f/]-)) 

+ ■■■1. (3.13) 



Here, we have only listed those terms that will actually contribute to the processes under study, 
when the mass-shell conditions for the momenta and polarizations of the photons are taken 
into account. We have however kept the numbering of the low-energy constants introduced in 
|27|, but we have, for convenience, changed their normalization by an overall factor l/Air"^. The 



notation is otherwise as in , except for the fact that the source x does not contain the factor 
25, see Eq. (0). 

The last piece we need for a full one-loop computation of the amphtudes ( ^751) is £(^2)- 
It counts as order 0{p^) in the standard case, and is not available from the existing literature. 
These contributions, which are order Oijn?) corrections to £(^o)' expected to be small 
in the two-flavour chiral expansion, and will be parametrized appropriately in the one-loop 
expressions of the amplitudes A^ and A^ given in the next section. The determination of the 
combinations of low-energy constants that enter these amplitudes will be discussed in Section 
5 below. 

When studying a given process one also needs to take into account contributions from 
pion loops, which produce divergences that are eliminated by a renormalization of the low- 
energy constants of the effective lagrangian. We have computed this divergent part of the 



*For a review, see 
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one-loop generating functional in the even intrinsic parity sector using standard heat-kernel 
techniques, and we have then performed the corresponding renormalization of the low-energy 
constants in the same dimensional renormalization scheme as described in |jT5|, 0. Thus the 
low-energy constants display a logarithmic scale dependence denotes generically any of 

these renormalized low-energy constants) 

X(/x) = + ■ ln(/x7/x) . (3.14) 

The full list of the resulting /3-function coefficients Tx is given in Table 1 below.Q 



X 


F'-Tx 


X 


F^-Tx 


X 


Tx 


A 






AB 


h 


1 

3 


ZP 




('^) 
Pi 


-AB{A + Z^) 


k 


2 
3 


ho 





(2) 
P2 


-AB{A - 3Z^) 


h 


1 

6 


h 




(2) 

Ps 


2B{A + 3Z^) 


k 


1 

3 






(2) 
Pi 


2S(3A + Z^) 




1 

12 






Ps 


AB{A - 2ZP) 







Table la: Scale dependence of the low energy constants of C^'^\ C^^^ and of C 



X 


F'-Tx 


X 




ai 


~2ZP 


ei 


-4^2 - 22(Z^)2 - 20AZ^ 


02 


-VIZ^ 




-4^2 - 16(Z^)2 - VIAZP 


03 







-12A^ - 6A{Z^y - ASAZ^ 


hi 


6(A + Z^) 


/i 


3^2 + 15(Z^)2 + UAZP 


b2 


-2(A + Z^) 


/2 


2^2 + 8(Z^)2 + lOAZ^ 


h 





/3 


AA^ + 24(Z-P)2 + 20AZP 


Cl 


2(A + 2Z-P) 


/4 


-6{ZP + A)ZP 


C2 





/5 


2^2 + 8(Z^)2 + lOAZ^ 


C3 





/I4 


AA^ + 8(Z^)2 + SAZ^ 


C4 


2A 


/l5 


-AA^ - 32{Zy - 28AZ^ 


C5 





/l6 


AA^ + 6(Z^)2 + SAZ^ 


/is 





/l7 


-iA{zy 



Table lb: Scale dependence of the low energy constants of £(^2,2) O'^d o/£(o,4). 

Notice that at order 0{p^), the low-energy constants of £'^2) and £(3) ^Iso need to be renor- 
malized. The corresponding counterterms, however, are of order (9(i?2) and 0{B), respectively, 

^ These results have also been established independently by L. Girlanda, private communication to M.K. 
and The renormalization of the Cf^ counterterms has, of course, already been obtained before in JlSt . 
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and they are gathered in the three last terms of Eq. ( p.8|) : in GChPT, renormahsation proceeds 
order by order in the expansion in powers of B. Ahernatively, one may think of Eqs. ( |1.6| ) and 
( p.7| ) as standing for the combinations C^"^^ + B'^C^^ 2) + BC^^ i) + BC'f^ 3^, respectively, 

with the corresponding low-energy constants representing the renormalized, scale dependent, 
quantities. We shall adopt the latter point of view in the sequel. 

We have not worked out the general structure of the divergent part of the one-loop 
generating functional in the odd intrinsic parity sector. For SChPT, this has been done in 



Refs. ra, 30, 34 



4 The One— Loop Amplitudes 

Having constructed the effective lagrangian in the preceding section, the computation of the 
amplitudes and at next-to-leading order is a straightforward exercise. We begin with 
the one-loop expression of the neutral amplitude A^ , whose structure at that level of the chiral 
expansion is given by ( p.8| ), with 



Ai {k,k';pi,p2,P3) 



and 



s-Ml 



-(/?-!) + -mt(/? - 3) - aM'J - AmH" - 700 



47r2F3l2 
+ -l(Ai + 2A2)(s-3M^)) 

TT 

+ Y2J{P12M2' Ml + \{a-l)Ml] 

+ ^[Ripl2, k ■ P12) + R{P12. k' ■ Pi2)M2 -M'^ + ha- 1)M2] 
+ ^[Ripls, k ■ Pis) + Ripls, k' ■ pr,)][p\z - Ml + \{a - l)Ml] 
+ YMpI,. k ■ P23) + R{pIs, k' ■ P23)]\pIs -Ml + ^{a- l)Ml] 



(4.1) 



A^{k,k';p^,p„p3) = ^^^R{pl„k-pu)[pl2 -Ml + ^{a- l)Ml] . (4.2) 
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This second amplitude, which was absent at tree level, is entirely generated by the pion loops. 
The numerator of the contribution which develops a pole at s = in the expression ( |4.1| ) of 
Ai is given by the product of A'""^'^'^, which is related to the on-shell amplitude of the 7r° — > 77 
decay through ^(vr" 77) = —ti{k, e, k', e')A'^°^'^'^ , times the amplitude ^°°'°°(p?2?Pi3 5P23) of 
virtual vr" — 7r° scattering (pfg + P13 + P23 = s + 3M^). The expressions, at order 0{p^) and 
order O(p^), respectively, of these amplitudes read 



A-" 



+77 



[l + mt + M;t' + mH" + 2m^a3 



(4.3) 



and 



aM. 



- + 4i(^i + 2A2)[(p?o-2M: 



2\2 



F2 ' F4^---^ ' -■^^LV/'12 

TV TT 



(pf3-2M,^)^ + (p^3-2M, 



r2\2i 



+ ^Jipl2 



1 

TT 

1 



+ ^Jipl3) 



(P 



23 



(4.4) 



The various parameters a, /5, Ai,2 and 700 involved in the expressions ( [4.1| ) and ( [4.4| ) are 
given in terms of combinations of the low-energy constants of the effective lagrangian and 
of chiral logarithms due to the pion loops. They read 



F2 



2mB + IQm^A 

+m^(81pS^^ + p?^ + 164pf + 2pf ) - 4M^m^(2) 

+ 16m^(16ei + 62 + 32/i + 34/2 + 2/3 + 72/4 + 603^) 



-8M^m'^{2bi - 262 
1 



as - 4ci; 



M2 



16772^2 



[2M^ + 102m' A + 264mM'] In 



12- 



1 A/f^ 



(4.5) 



= l + 2me(')-4m2(e(2))2 



2m2(3a2 + 2a3 + 46i + 262 + 4ci) 



48.2^n^^f + + 



(4.6) 
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487r2 ^ ' 



(4.7) 



7oo 



m\3a2 + 2a3 + 6b2 + I2ci) 



Ml 



32772^2 



M2 

ln^ + l)(a-l) 



(4.8) 



The parameters t, t' and t" contain the contributions from C^^^^y '^(eo) '^(42)' respectively. 
In particular, from the formulae (|3.12| ) and (|3.13|) we derive 



t 
t' 



32 , 
3 

-^(^2-2^3) 



(4.9) 



Finally, J{s) denotes the Chew-Mandelstam function pS, the usual scalar two-point loop 
integral subtracted at s = (for its expression, see Ref. fl^), whereas RijP'.k ■ p), which is 



related to the three-point scalar loop function, is given, for A;2 = g, by 



i?(p2^ k.p) = C{p^, k-p)- 



{k-pf 
4{k ■ p) 



[J{{k-pf)-J{p')] + 



327r2 



(4.10) 



with (a = Y^l - 4M2/p2, a' = AM^/{k - p^ 
167r'^C{p'^, k ■ p) 



a +1' 



A{k -p) 



ln^(^)] 



(4.11) 



For the charged amplitude , the general structure is more involved, see Eq. (|2.6| ), and 
at one loop we obtain 



/1^°^77 /100;+-/ 2 2 2 \ 

^f(A;,A:';po,P+,P-) = ^^+-'^+0'^-o^ 



s-Ml 

+ 4^Ff {^(^ + + r'^^ - - - 
+7+- + 37;. 



TT 7r 



15 



1 



f— ,J(pi_)[3pi_ + 4(a-l)M,1 



6F2 
1 



^V-o)[5/-o-14M^'-2(a-l)M^] 



(4.12) 



-^J(p%)[4M^ - i^J(p2_^)[4M^ 

TT TT 

-^J((fc'-Po)^)[4M,'-(A:'-Po)l}, 



(4.13) 



^3 {k,k';po,p+,P-) 



-^^R(pl_-k.p^_)\pl_ + ^{a-l)Ml] , 



(4.14) 



A {k, k';po,p+,p_) 



P+o 



/>+o + 2M^ + i-4^(4M^^-p^o) 



k-p+o 



{AM^-7p% + 6k-p+o) + 2 



2k -p^ 
k ■ P+ pIo 
k-p+ok- p+o 



(p% - iM^) 



+R{pI,, k ■ p+o)[ - + 2A; • pAt^ " 1) - - 1)^^] 



P+o' 



(4.15) 



A {k,k';po,P+,P-) 



247r2F5 



J((A;'-p_)2)[4M^(A;'-p_)V^>+o)(4M^p%)- 



247r2 



^■P4 

(4.16) 



A (/c,/c';po,P+,P-) = 



47r2F5 1 6 
+4 



P%-4M2-6A;-p+ 



A: - p+o 



(/^•P+o)^ 



-J ■ P+o 
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+ 



1 k ■ p+ 
327r2 k ■ p+o 



{k'-p_y + 



967r2 V3 



-k ■p+ -p 



+0 



(4.17) 



The amplitude again contains a contribution with a pole at s = M^, which is given 
by the product of the O(p^) vr^ 77 amplitude ^'^"^tt times the (off-shell) ti^tt'^ tt+tt^ 
amphtude A'^^'~^~ , with 



00;- 



-/ 2 2 2 \ 

(P+-,P+o>P-o) 



/5 .^2 



F2 
1-^1/2 

TT 

1 , 



3 

2M!)2 



3F2 

+ ||[fc-2il/a^ + (pi 



6F4 
1 



+ 



12F4 

TT 

1 



12Fi 



ApI-) 
ApIo. 

Ap'-o. 



4(P^ 



2Mlf] 
4 



3(p^o - ^M,' - l^^Mlf + (p^ 



3(P% - \mI - \aMlf 



where + p% + p% = s + 3M2. 

The remaining parameters 7+_ and 7;^^ which appear in the amplitudes ,4f , 
contain the contributions from the low-energy constants and chiral logarithms : 



-pWpl^-mi 

^{P\--P\,){P\-^MI 

(4.18 



7+- 



-m2(a2 + 262 + 4ci) 



M2 M2 



-) +m2(57 



+- ) 



(4.19) 



7;_ = 8m(2Ai2 - A13) + 8M2(A7 - As) + em^as 

32772^2 /i2 ^ 



(4.20) 



With the expressions given above, it is straightforward to check that the three conditions 
listed before Eq. ( |2.9| ) as well as Eq. ( p. 10 ) are satisfied, which provides a non-trivial check of 
our calculation. The isospin relation ( p. 9]) implies that the condition 



2 + 37+_ + 7oo + 



Mi 



w?{2a^ - 3t") 



(4.21) 
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must hold. This requires that the contributions of the -^^^2) counterterms to 7^ , which we 

have denoted as m^57^ in the expression ( [4.11j ), have to satisfy 



3t" + 3(574 







(4.22) 



The contributions of the 2) counterterms to 



_, m^57^_ in the expression ( |4.20| ), are not 
constrained by isospin symmetry. 

Upon using the information provided by Table 1, it is straightforward to check that a, (3 
and 7oo are scale independent by themselves (the parameters Ai and A2 were directly expressed 
in terms of the scale independent quantities h and I2 defined in Ref. [|l5l). In order that the 



amplitudes A'"°^'^'^, A'^ and A'~^ be independent of the subtraction scale /i, the parameters t, 
t', t", 7+_ and 7^_ must be separately scale independent. This requires that 57+_ and (57^_ 
themselves are scale independent. Since we have not worked out the structure of the one-loop 
divergences of the generalized generating functional in the odd intrinsic parity sector, we could 
not perform these checks explicitly. 

^From the above formulae, one may infer the expressions of the amplitudes in the stan- 
dard case |2^, Since only the result for the neutral amplitude A^ was displayed explicitly 
in Ref. ||2^, and the expressions of the amplitudes Af in the charged channel are only available 



from the unpublished work p6[, we describe in some detail the necessary steps to obtain them. 
Their general structure is of course unchanged, the differences occur only in the expressions of 
the various combinations of low-energy constants that are involved. In particular, the contri- 
butions from >C(o,3), '^.{2,2), -^(0,4)5 and >C(^2) relegated to higher orders. For the remaining 
constants, the correspondance with the usual SChPT notation is given as follows 



. ^(2) 



/5st 
7oo,st 

7+-,st 



167r2F2 



M2 

(^4 + In 



1 + 
, 



327r2F2 
M2 



(4/4 - 3/3 - 1) 



87r2F2 



{k - 1) , 



Ml 



727r2F2 



7V- 



,st 



2Ai2- A 



13 ^ 



Ml 



2B 



In 



Ml 



327r2F2 /i- 



(4.23) 



We have checked that upon substituting these expressions into the one-loop amplitudes A'^ 
and A*" given above, we recover the results of the standard case, up to the contributions from 



the counterterms contained in t and t', which were not included in Refs. [E^, EG 



18 



5 Counterterm Estimates and Numerical Results 



In order to make numerical estimates for the cross sections based on the O(p^) amphtudes, we 
first need to fix or estimate the values of the various counterterms involved. 

i) Ai, A2 : At order 0(p''), these parameters are related to h and I2 through Eq. (|4.7| ). The 
values of these low-energy constants in the standard case have been the subject of numerous 
studies in the hterature [|T5|, |3^, ^ A first determination at order 0{p'^) was given 
in Ref. [0, using information from the D-wave vr-vr scattering lengths. The corresponding 



values, taken from a recent numerical re-analysis [^], are 

[i,GL = -2.15 ±4.30 , /2,GL = 5.84 ± 1.72 , (5.1) 

leading to 

Ai^GL = (-7.35 ± 9.06) X 10"^ , Aa.cL = (10.57 ± 3.63) x 10"^ . (5.2) 
The parameters Ai and A2 can also be determined directly, via a set of rapidly convergent sum- 



rules ^T], from the knowledge, at order two loops, of the vr-vr scattering amplitude A{s\t,u) in 
GChPT |12|, and from medium energy data on vr-vr phase shifts. The values obtained this way 
correspond to a determination at order 0{p^). They depend only very weakly on the values of 
the parameters a and jS when the latter are varied within the ranges specified below, and read 

Ai = (-6.1 ±2.2) X 10"^ , A2 = (9.6±0.5) X 10"^ . (5.3) 



These values are compatible with those given in eq. (|5.2|) , but are affected by much smaller 
error bars. The analysis may even be refined in the standard case, using the information on 
the SChPT two-loop vr-vr amplitude obtained in Ref. [^], leading to the following values ^9 



Ai,st = (-5.7±2.2) X 10-^ A2,st = (9.3±0.5) X 10-^ (5.4) 

ii) a, P : At leading order, a is directly correlated to the size of the condensate, see Eq. (|1.12| ). 
As such, the value of a is not predicted by GChPT, but remains a free parameter, that can a 
priori be varied in the range 1 ^ a ^ 4. At order (9(p^), the relationship between a and the 
ratio 2mB/Ml becomes more complicated, as shown in Eq. ( [4.5|) . The corresponding next- 
to-leading and next-to-next-to-leading corrections have been estimated in Ref. [0 (see in 
particular the figures 8 and 10 in that reference). Notice that the analysis of Ref. |T2| was 
done within the framework of 5'?7(3)l x S'f/(3)R chiral perturbation theory. Working with 
only two light flavours as in the present paper might further reduce the uncertainties in the 
correspondance between the value of a and the size of the condensate, due to the lower number 
of unknown counterterms involved, and due to the absence of large contributions from the chiral 
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logarithms induced by the kaon loops ||3^. However, the results of Ref. are sufficient for our 



present purposes, and we shall not pursue that matter further. Once a is given, the parameter 
13 is also constrained by low-energy tt-tt data. The correlation between a and (3, which results 



from the Morgan-Shaw universal curve has also been studied beyond leading order in 
I2| , and is summarized in Fig. 6 of that reference. For the subsequent numerical analyses, 
we shall take the values given in Table 2 below. The values given in the second line of this 
table correspond to the one-loop values a^t and j3st of the standard case, which follow from 



the expressions given in Eq. ( |4.23| ), and from the values = 2.9 ± 2.4 [0, U = 4:-4: ± 0.3 
Two points are worth being remembered. The first is that, independently of the value of a, /3 
stays close to unity. The second point is that the standard case allows to make a very precise 
prediction for the value of a at order 0{p^), viz. agt = 1-06 ± 0.06. Furthermore, this value 
is barely affected by the corrections at order O(p^) : The analysis of Ref. |3^, based on the 



results of [^, gives ctst = 1.07±0.01 and Pst = 1.105±0.015 at next-to-next-to-leading order. 



Therefore, any significant deviation of the value of a from unity would provide evidence for a 
departure from the standard scenario of chiral symmetry breaking with a strong condensate. 

iii) t, t', t" : The constants t, t' and t" appear in the expression of the 7r° — 77 amplitude 
( [4.3[ ). The uncertainty on the experimental value of the decay rate r(7r° — » 77) = 7.74±0.56 
eV, only yields a very weak constraint on the combination that appears in A^'^"*'^'^' , viz. mt + 
M^t' + m'^t" = (0.0 ±3.6) X 10^^. This is comparable to the estimate one would obtain through 
naive dimensional analysis MM. In addition, as shown in Ref. [EH], isospin breaking effects can 



,0 



be sizeable in A" , Further information may be obtained by making use of the sum-rules 



considered in Ref. The corrections due to t' were however not taken into account there. 



but the analysis of is easily extended to the more general situation. In the generalized 
case, a similar set of sum-rules can be established, but they do not yield complete information 
on the three constants t, t' and t" . We have summarized this analysis in the Appendix for the 
interested reader. Here, we only quote the values that we shall use in the sequel (the estimate 
for t" follows from naive dimensional analysis), 

mt = (6 ± 12) X 10"^ , Mlt' = (-3 ± 3) x 10"^ , mH" ~ ±1 x 10"^ . (5.5) 

iv) 7oo, 7h , 7^_ : The main difficulty in obtaining numerical estimates for these parameters 

comes from the lack of knowledge, in the generalized case, on the contributions from the low- 
energy constants of £^ 2)- spirit of Vector Meson Dominance (VMD), we have estimated 
the contributions of vector mesons to the two amplitudes and A'" , as done in Ref. [25 
for the standard case. This way, we find that t, t' and t" receive no contribution, whereas the 
contribution to 700, 7+- and 7^_ read 

3 

7oo|vMD = 7+-|vMD = , 7+_|vMD = ~^J^ 0-024 . (5.6) 
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We take (|5.6|) as the values of these constants at the scale /z ~ My=770 MeV, and estimate the 
error associated to the VMD approximation and to the presence of the low-energy constants 
from £^ 2) ill the expressions of 700 and of 7+^. by varying the scale /i of the corresponding 
chiral logarithms between 500 MeV and 1 GeV. The resulting values for 700 and 74— are shown 
in Table 2. In the case of 7+_, we obtain a constant value which, to a very good precision, is 
compatible with zero. 



a 


/? 


700 X 10^ 


7+_ X 10^ 


1.06 ±0.06 


1.103 ±0.008 





-3 


1.5 


1.06 ±0.06 


8±3 


-5.7 ±0.8 


2 


1.07 ±0.06 


16 ±5 


-8.4 ± 1.6 


2.5 


1.08 ±0.06 


24 ±8 


-11.0 ±2.4 


3 


1.11 ±0.06 


32 ± 10 


-13.7 ±3.2 



Table 2:Values of (3, 700 and 7+_ for different values of a. 



With the above inputs at hand, we may now consider a few numerical applications. In 
Fig. 2, we have plotted the tree-level and one-loop cross sections for the charged channel, 
(T^gg(s,a) and (T'"(s,a), in the threshold region, 3Mt, < ^/s < 0.5 GeV, where we expect the 
one-loop expression of the amplitude to be reliable, and for different values of a. We have used 
M^-i = M^o = 135 MeV in the amplitude and the experimental values, as quoted in P|, in the 
phase space integrals. 

Let us first concentrate on the standard case (a ~ 1) which has been discussed before 
in the literature. In the second half of the energy region that we have shown, the correction 
as compared to the tree level cross section amounts approximatively to a factor of two, if we 
consider the central values. This agrees with the unpublished result ||2^, but is much less 
than previously found in . On the other hand, the error bars induced by the uncertainties 
attached to the various conterterm contributions that enter the 0(p^) amplitude are important. 
A closer analysis (see also below) reveals that the main contribution comes from the uncertainty 
on the value of Ai given in (|5.4|) . This shows also that the sensitivity of the cross section on 
the 0{p^) counterterms becomes already sizeable even at such low energies. Coming now to 
the dependence with respect to a, we see that here also the higher order corrections have a 
deep influence and upset the situation that prevailed at tree level. For a > 2 the one-loop 
cross sections are suppressed as compared to their tree-level values, and become even smaller 
than o"'"(s, 1) as the energy increases. Unfortunately, the present theoretical error bars make it 
difficult to disentangle the different situations in practice from the knowledge of the total cross 
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section alone. 
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Figure 2 : The cross section a{s,a) (in logarithmic scale) for^'j tt+tt^tt^ as a function of 
the center of mass total energy, for three values of a. Also shown are the corresponding curves 
for the tree level cross section atree{s, o:). 



The cross sections a(^^^{s,a) and a^{s,a) in the neutral channel have been plotted in 
Fig. 3. Whereas cr^ee('^> ~ o;^) corrections are seen to have an even more drastic influence 
on the behaviour of the cross section as a function of energy than in the charged channel. 
Unfortunately, as far as the dependence on a is concerned, the picture is again totally blurred 
by the uncertainties, which, for the sake of clarity, we have not shown, but which are even more 
important than in the charged case. 

In both channels, the origin of the large error bars is a consequence of the highly de- 
structive interferences between the various amplitudes. In the charged case, this interference 
was already present at tree level, and is even accentuated by the loop effects. In the neutral 
case, the strong a dependence of the single amplitude that contributes at lowest order is, in 
a similar way, washed out by the interferences between the two one-loop amplitudes Ai and 
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A2 ■ In order to illustrate this point, we show, in Fig. 4, how the neutral cross section looks 
like if only the amplitude Ai is considered. The importance of the interference effects with the 
second amplitude A2 appears clearly upon comparing Figs. 3 and 4 (for the standard case, a 
similar observation was already made by the authors of Ref. |^). Unfortunately, we have not 
found a simple way to extract the contribution from A^ only : for photons with perpendicular 
polarizations, the two amplitudes contribute, and the destructive interference between them is 
again at work, while for photons with parallel polarizations, only the very small contribution 
from A2 , which vanishes at tree level, is singled out. We have tried to investigate whether 
looking at more refined observables allows to reach better perspectives from this point of view. 
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Figure 3 : The cross section a{s,a) (in logarithmic scale) for 77 — tt^tt^tt^ as a function of 
the center of mass total energy, for three values of a. Also shown are the corresponding curves 
for the tree-level cross section atree{s,a). 



We have, for instance, looked at the invariant mass distribution of the two charged pions in 
the 77 TT+TT^Tr" channel. The result is shown in Fig. 5. As one may observe, the error bars 
are much less important than for the total cross section, and the different values of a can be 
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distinguished over a substantial portion of the energy range that has been considered. Actually, 
analyses of this type usually require sufficiently high statistics, which also represents a problem 
in the present case. 
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Figure 4 : The cross section cr^]v(s,a) for — > 7r°7r°7r° obtained by taking into account the 
contribution from the amplitude alone, for different values of a. 

Indeed, in both channels, the cross sections are very small at low energies, orders of 
magnitude below, for instance, the corresponding cross sections for the 77 —>■ tttt processes. 
We have therefore also estimated the numbers of events that could be expected at an e"*" — e~ 
collider for the two-photon total invariant mass ^/s below a maximal energy E'max- As typical 
examples, we have considered two instances of symmetric e"*" — e~ colliders. The first case 
corresponds to the Daphne 0-factory with a total beam energy of i?boam = 510 MeV, 
and a nominal integrated luminosity of 5 x 10^ nb~^ per year. The second case concerns a 
r-Charm Factory configuration [^, with a beam energy four times as large as for Daphne, 
and a design integrated luminosity of 10^ nb~^ per year. The numbers of events are obtained 
upon convoluting the above cross sections with the corresponding photon luminosities quoted in 
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1 47 1 (since we are not interested in resonant rj production, whenever necessary we avoided the rj 
peak by applying a cut on m^y such that only the events with m^^ < — A or m^^ > + A, 
with A = 20 MeV are accepted), 



dL 



77 



dm. 



'77 



77 



' ^ In ^^^^'^ ^ ^ f 



)^[-(2 + ^^)^ln^-(l-^^)(3 + z^)]. 



(5.7) 



where m-y^ = is the photon-photon center of mass energy, z = m^^/2£'beam, and rrie is the 
electron mass. For the luminosity quoted above, the expected total number of 77 tt^tt^tt^ 
events per year at energies ^/s < 0.6 GeV is around 5 ± 1 for Daphne (independently of a) and 
even less in the neutral case. In the case of a r-Charm Factory, the total number of events 
becomes sizeable, and the results are shown, for the two modes, in the fourth column of Table 3, 
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Figure 5 : The histogram of the distribution of the invariant mass m,r+7r- of the two final 
charged pions in the process 77 tt+tt^tt^ for different values of a and ^/s < E"^"-^ = 0.6 GeV. 
Also shown are the curves correponding to the same distribution, but from the lowest order 
amplitude alone. 
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Table 3: Number of events for 77 — > tt'^tt vr*^ and /or 77 Tr'^vr'^Tr^ at a r-Charm Factory as 
a function of the maximal energy E^^^GeV) , and the principal sources of error. 

for different choices of -Emax (second column) and for different values of a (third column). Also 
shown are the corresponding uncertainties A(7^events) (fifth column), while the remaining 
entries show various sources of contributions to the total error A (^events). The latter was ob- 
tained upon adding the individual contributions in quadrature. As mentioned above, the main 
source of error comes from the uncertainty on the value of Ai. A sizeable but not necessarily 
drastic reduction of the latter would already allow to distinguish the standard case of a strong 
condensate (a ~ 1) from situations where spontaneous breakdown of chiral symmetry would 
be triggered by a much weaker condensate {a ^ 2). 



6 Summary and Conclusions 

In the present paper, the amplitudes of the processes 77 — > 7r°7r°7r° and 77 7r''"7r~7r° have 
been computed in the framework of SU{2)i^ x SU{2)^ Generalized Chiral Perturbation Theory 
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to 0{p^) precision. The corresponding generating functional has been constructed exphcitly in 
Section 3, and the structure of its divergences in the sector of even intrinsic parity has been 
analysed. The resulting amplitudes, worked out in Section 4, satisfy the isospin relations that 
we have established in Section 2 (to the best of our knowledge, these relations have not been 
discussed previously in the literature). When restricted to the standard case, specified by the 
choice of parameters as indicated in Eq. (|4.23| ), we recover the results obtained by previous 



authors |]25|, both in the neutral and in the charged case. Finally, we have estimated 
the counterterms that enter the one-loop amplitudes and we have performed some numerical 
analyses in Section 5. We have, in particular, shown that the error bars associated to the 
counterterm estimates become important, especially in the neutral channel, as a consequence 
of highly destructive interference effects between the various amplitudes that build up the total 
cross sections. 

We have also considered the possible detection of these processes at Daphne and at a 
r-Charm Factory. Unfortunately, and precisely because of these large interference effects, the 
expected number of events is rather discouraging in the first case. Depending on the actual 
values of the counterterms and on a, it is hard to expect more than ^ 5 events per year with 
total invariant mass lower than 500 MeV. The number of events increases substantially when 
allowing larger invariant masses, but at the expense of working in an energy region where the 
0{p^) expressions are probably less reliable, since higher order terms can become important. 
The computation at 0{p^) would allow a better control of the cross sections at larger momenta, 
thus allowing a substantial increase of the number of events already at Daphne. The required 
amount of work seems excessive, though, and would make sense only if conducted in parallel 
with a better determination of li (for instance, from a two-loop analysis of Ki^ decays), by 
far the main source of theoretical uncertainties at present. We rather expect interesting and 
realistic prospects in this field to come from future machines, like the r-Charm Factory, which 
run at higher energy and higher luminosity. 
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Appendix 

In this Appendix, we give a brief description of our analysis of the anomalous counterterms 
A2, A^, A4 and Aq which is based on the approach of Ref. [^. The starting point are the 
invariant amplitudes nyyp(p^, g^, r^) and Ilyyp{p'^,q'^,r^) of the vector-vector-pseudoscalar 
three-point correlation functions in the three flavour chiral limit (we take the definition given 
by Eqs. (3) and (4) of |^). At = = 0, the loop contribution vanish, and one has (notice 



the absence of the pion pole in the second equality) 

nv.yp(0, 0, r^) = + t\[8A4 - 16So(A2 - 2A^)] , 

nV(0,0,r2) = -L[8A4 + 24A6] . (A.l) 

In the chiral limit, the counterterms from £^ 2) do not contribute, so that the above result holds 
both in SChPT and in GChPT. Although the contributions from the low-energy constants A2 
and ^3 from were omitted in |^, one may follow the same steps as described there (we 
have also kept the same notation), and end up with the following set of sum-rules : 

^[8A-165o(^-2A3)] = -A(4 + ^f^ftan9'^(/-^^n , (A.2) 



and 



[24A + 16Bo{A2 - 2As)] = --^V6A{v' - 77) • (A.3) 

In addition, one needs to know the expression of the amplitude of the two-photon decay of the 
r], A{'ri — i> 77) = —ti{k, e, k', e')A^^^^, which, in analogy to the 7r° 77 amplitude A'""^'^^ of 
Eq. (f4.3|) , may be written as 

A"^^^^ = [— + -m(5 - 2r)t + — m(l - r)A6 + MH' + m^t"{r) + -m'il + 2r'')a3\ . 

AV^tt'^Ft, Fri 3 3 3 

(A.4) 



In this last expression, t and t' are related to A2, A^ and A^ according to Eq. ( |4.9| ), while r 
stands for the quark mass ratio mg/m, and m'^t"{r) denote the corrections coming from >C^2)) 
which can be of order (9(m^), 0{mms), and 0{m1). 

In the standard case, the two last terms in Eq. ([A.4| ) would appear only at higher orders, 



and one may thus proceed as described in Ref. At the order we are working, the quark 
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mass ratio is then given as r^t = r2 = 2M|-/M^ — 1 ~ 25.9 |16|. Using the numerical values 



given in one obtains, from Eqs. ( |A.2| ) and ( [A .31 ), respectively, 

rnA^,st - Ml{A2 - 2A^),, = (-5.9 ± 1.8) x 10"^ (A.5) 

and 

3mA,st + Ml{A2 - 2A3)st = (-2.1 ± 0.4) x IQ-^ . (A.6) 
Upon using the experimental rate for — >■ 77 P], we then determine the combination 

MliA^ - 2A3)st = (1 ± 1) X 10-3. (A.7) 

Adding errors in quadrature, the previous results give 

mt,t + = (1.7 ± 8.2) X IQ-^ , (A.8) 

which is four times more accurate than the value obtained directly from the experimental rate 
of 7r° — i> 77. For the separate pieces, we obtain 

mUt = (4.4 ±10.8) X 10-^ M^t = (-2.7 ±2.7) x 10"^ . (A.9) 

In the generalized case, the analysis may, unfortunately, not be pursued quite that far. 
The main drawback are the corrections from £^ 2) 5 which in particular produce potentially 
large (9(m^) corrections to the ^ 77 decay amplitude, and on which the sum-rules (|A.2| ) 
and ( |A.3| ) give no information. If one restricts the analysis to the order 0{p^) precision, then 
the contributions from C'^qq^ are also absent, and the situation becomes even simpler than in 
the standard case, since (c/. Eq. (|1.4|) ) the left-hand sides of the sum-rules (|A.2| ) and (|A.3| ) 



now only involve ^4 and Aq, respectively. Keeping in mind that r is now a free parameter 



(a and r are however related, see 0]), and taking the necessary inputs from the 0{p ) 
determination of ^4 and Aq reads 

_ iVc f \X{r)Ml-il-X{r))'M'J \'^ Ajn' ^ jj) , 
""^^-"32 1, M2-M| ) M;^(vr^77) ' ^ ^ 



where the ellipses stand for higher order corrections, A(r) = 2(r2 — r)/(r^ — 1) and Aqmo = 
(3M^ — 4M|; + M^)/M^ ~ —3.6. For r = r2, the above expressions give values compatible with 
the previous SChPT analysis. Furthermore, as r decreases {i.e as a increases), the variation 
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of frit stays within the bounds given in Eq. ( A. 9 ). On the other hand, in the extreme case of 
a vanishing condensate, the £^ q-j contributions to the two sum-rules also disappear, and the 



expressions ([A .101 ) and ( |A.11| ) become exact at order 0{p^). Finally, if we use the two sum- 



rules in order to express and Aq in terms of A2 — 2A3 in the expression (|A.4|) , we obtain 
an estimate of a combination of A2 — A3, m?t"{r) and rfi^a^, which is not very sensitive to the 
value of r and compatible with the value (|A.7|) obtained in the standard case for M^(A2 — 2^43). 



Thus, within reasonable error bars, the values of t and t' can be taken independent of r. For 
the numerical analyses presented in the text, we have used 

mt = (6 ± 12) X 10^^ , Mlt' = (-3 ± 3) x 10"^ . (A.12) 
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